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Abstract 

We provide an explicit structure of the charged lepton mass matrix which is 2-3 
symmetric except for a single breaking of this symmetry by the muon mass. We 
identify a flavor sjTumetric limit for the mass matrices where the first generation 
is decoupled from the other two in the charged lepton sector while in the neutrino 
sector the third generation is decoupled from the first two generations. The leptonic 
mixing in the symmetric limit can be, among other structures, the bi-maximal (BM) 
or the tri-bimaximal (TBM) mixing. Symmetry breaking effects are included both 
in the charged lepton and the neutrino sector to produce corrections to the leptonic 
mixing and explain the recent ^13 measurements. A model that extends the SM 
by three right handed neutrinos, an extra Higgs doublet, and two singlet scalars is 
introduced to generate the leptonic mixing. 
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1 Introduction 



We now know that neutrinos have masses and just hke the quark mixing matrix 
there is a leptonic mixing matrix. This fact has been firmly estabhshed through 
a variety of solar, atmospheric, and terrestrial neutrino oscillation experiments [1]. 
We parametrize the neutrino mixing matrix, Upmns, as follows [2]: 

(C12C13 S12C13 si3e~'^ \ 

-S12C23 - Ci2S23Si3e*'^ C12C23 - Si2S23Sl3e*^ S23C13 K, (1) 
S12S23 - Ci2C23Sl3e*'' -C12S23 - Sl2C23Sl3e*'' C23C13 / 

where S13 = sin ^13, C13 = cos ^13 with ^13 being the reactor angle, S12 = sin ^12, 
C12 = cos 6112 with 9i2 being the solar angle, S23 = sin 6*23, C23 = cos ^23 with ^23 being 
the atmospheric angle, 6 is the Dirac CP violating phase, and K = diag{l, e^"^^, e^'^'^) 
contains additional (Majorana) CP violating phases (pi, 4>2- We ignore the Majorana 
CP violating phases in this work. 

Unlike the CKM matrix that can be thought of as a perturbation about the 
identity matrix the leading term in the leptonic mixing contains large mixing angles. 
Some examples of the leading order mixing matrix are the bi-maximal mixing [3] 
and the tri-bimaximal mixing [4]. However, current experiments indicate deviations 
from these standard zeroth order forms. 

For instance, recent results from the T2K [5] and MINOS [6, 7] experiments 
have indicated a large reactor angle ^13 for neutrino mixing. At the 90% C.L., 
T2K gives 0.03(0.04) < sin2 2^i3 < 0.28(0.34), with zero Dirac CP phase, Sd, for 
normal (inverted) hierarchy. The MINOS group gives 0.01(0.026) < sin^2^^i3 < 
0.088(0.150). There are already several papers that have attempted to explain the 
resent 6*13 results [8, 9]. Large 6'i3 was anticipated in Ref. [10]. 

The leptonic mixing arises from the overlap of matrices that diagonalize the 
charged lepton and the neutrino mass matrices. Many approaches to studying the 
leptonic mixing start in the basis where the charged lepton mass is diagonal. Our 
approach to obtaining the leading order leptonic mixing as well as deviations from 
it starts from the charged lepton sector. A recent attempt to understand Oi^ from 
the charged lepton sector can be found in Ref. [8] and in the past corrections to 
the leptonic mixing from the charged lepton sector were considered in Ref. [1 1] . An 
approach to suppress flavor changing neutral current effects (FCNC) in the quark 
sector, based on shared flavor symmetry, was proposed in Ref. [12]. As an example 
of this shared symmetry the decoupled 2 — 3 symmetry was used for the down 
quark sector to suppress FCNC effects and explain anomalies [13] observed in the 
B meson system. In the decoupled limit the first generation is decoupled from the 
other two generations. We extend this decoupled 2 — 3 symmetry to the charged 
lepton sector. This is a reasonable extension given the fact that the down quark and 
charged leptons exhibit similar hierarchical structure and they may be combined in 
representations of GUT groups. 



1 



One of the central ideas of this approach is the requirement that the mass ma- 
trices, in a symmetric hmit, be diagonahzed by unitary matrices composed of pure 
numbers independent of the parameters of the mass matrices. This is similar to 
the idea of form diagonahzable matrices discussed in Ref. [14]. If one starts with 
a 2 — 3 symmetric mass matrix for the charged Icpton sector and requires it to be 
diagonahzed by unitary matrices of pure numbers one recovers the decoupled 2 — 3 
symmetry. In the neutrino sector we assume the third generation to be decoupled 
from the first two generations. With real entries in the neutrino mass matrix it is 
diagonahzed by a rotation matrix and the resulting leptonic mixing has a — r 
symmetry. Requiring the mass matrix to be diagonahzed by pure numbers can lead 
to, among other structures, the BM and the TBM leptonic mixing. 

To generate the mixing matrices in the charged lepton and the neutrino sector, 
we present a Lagrangian that extends the SM by three right handed neutrinos, an 
additional Higgs doublet and two singlet scalar fields.^ The Lagrangian uses the 
same class of symmetries as has been used in Ref [16]. However, the structure as 
well as the phenomenology of our model is very different from the above mentioned 
papers. The Lagrangian is constructed to have a 2 — 3 symmetry, along with 
two additional Z2 symmetries Z| and ■ The neutrino masses and mixing are 
generated through the usual see-saw mechanism. The presence of the x Z2 
symmetries lead to the decoupled 2 — 3 symmetry in the charged lepton sector and 
fixes the interactions of the right handed neutrinos with the singlet scalar fields. 
The presence of the Z2 symmetry forces the neutrinos to acquire Dirac masses by 
coupling to a second Higgs doublet which has a different Z2 transformation than the 
usual SM Higgs doublet that give masses to the charged leptons. The full Lagrangian 
is symmetric under the product of the Z2 symmetries, x Z| x Zl'. 

The neutrino masses and mixing arise when the Higgs doublets and the singlet 
scalars acquire v.c.v's and break the symmetries of the Lagrangian. The leptonic 
mixing is predicted to be of the bi-maximal type when both the singlet scalars 
acquire the same v.e.v. If the v.e.v of the second Higgs doublet is small enough 
~ MeV then the see-saw scale as well as the masses of the singlet scalars can 
be in the TeV range. To obtain the TBM mixing one has to use different fiavor 
symmetries. 

Symmetry breaking is introduced in the charged lepton sector by higher dimen- 
sional operators that break the decoupled 2 — 3 symmetry but generate a 2 — 3 
symmetric mass matrix except for a single breaking generated by the muon mass. 
In the neutrino sector, symmetry breaking is introduced by breaking the alignment 
of the v.e.v's of the singlet scalars by terms in the effective potential. The correc- 
tions to leptonic mixing go as ~ ^ where v is the v.e.v of the SM Higgs and a; the 
scale of the singlet scalar v.e.v's. If a; ~ TeV then the corrections to the leptonic 
mixing are enough to explain the experimental observations. 

The paper is organized in the following manner: We begin in Sec. 2 with a 

^Recent motivations for considering two Higgs doublet models can be found in Ref. [15] . 



2 



discussion of the flavor symmetric limit that leads to among other structures the 
BM and TBM mixing. In Sec. 3 we present the Lagrangian to generate the mixing 
matrices in the symmetric limit. In Sec. 4 we study the effect of symmetry breaking 
in the charged lepton and neutrino sector to generate the realistic leptonic mixing 
matrix. In Sec. 5 we show the numerical results due to the symmetry breaking, and, 
finally, in Sec. 6 we conclude with a summary of the results reported in this work. 



2 The Leptonic Mixing in the Symmetric Limit 

We start with the charged lepton sector, and assume that the Yukawa matrix is 2 — 3 
symmetric [17]. The Yukawa couplings of the charged leptons are given by 

(hi 1x2 —^12 
h2 I22 I23 I . (2) 
—lu I23 I22 

The above Yukawa matrix can be diagonalized as 



(3) 

where the mixing angle 9 is determined by the positive solution to 

tan^ = '^ f^'' — ■ (4) 

I22 — I23 ~ ^11 i v {I22 — I23 ~ hiY + 8^12 
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The eigenvalues of are + I22 - ks ± V {hi - I22 + kzY + 8^12] and I22 + ^23- 
According to our assumption, the elements of the matrix that diagonalizes Y^ must 
be pure numbers in the symmetric limit. It is clear that we can achieve that by 
setting /12 = (6* = 0) in Eq. 4. This generates the decoupled 2-3 symmetry [12], 
as the flavor symmetry in the charged lepton sector in which the flrst generation is 
decoupled from the second and third generations. 

One can represent the Yukawa matrix with the decoupled 2-3 symmetry by Y^ 

as 

(hi \ 

Y^, = i/22 . (5) 

\ 2^23 2^22/ 

This Yukawa matrix 123 is diagonalized by the unitary matrix W^23 given by 

no \ 

< = 71 • (6) 
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Note that this matrix differs from the one in Eq. 3 in the hmit ^ = by an irrelevant 
diagonal phase matrix. Writing the diagonalized Yukawa matrix as Y.^^ we have 

'hi 

Y^sa = WgY^,Wl,^{ |(/22 -/23) |. (7) 



The charged lepton masses are given by 

J. ^1 ; 
rUe = ±—^lu, 

v2 

, vi {I22 - hs) 
= ±^(^. (8) 



Since << rur there has to be a fine tuned cancellation between I22 and I23 to 
produce the muon mass. Hence, it is more natural to consider the symmetry limit 
^22 = hs which leads to = 0. The Yukawa matrix which leads to the zero muon 
mass within the decoupled 2-3 symmetry is 

hi 

y^i = I \It I . (9) 

In the neutrino sector we assume that, in the symmetric limit, M.^, has the 
general structure 



= \ d h , (10) 





where all the parameters are real. This can be diagonalized by the matrix 

(C12 S12 0' 
S12 -C12 
1 

S12 = sin 6112, 

C12 = cos 6*12, (11) 



where 

(a -6) 



2d 

tan 2^12 = T-. (12) 
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We can then calculate Upj^^g as 



U p 



PMNS 



(13) 



with 




(14) 



Cl2 







) 



UpMNS = I "73^12 75C12 



1 



1 



(15) 



which is just the ji — t symmetric leptonic mixing. If we require 612 in Eq. 12 to be 
independent of the parameters a, h and rf, then, we either have a = h which leads to 
O12 = 7r/4 and generates the BM mixing or d = k{a — h) and in particular we obtain 
the tri-bimaximal mixing with k — -\/2. Hence, by choosing a — b, the neutrino 
mass matrix is given as 



3 The Lagrangian in the Symmetric Limit 

In this section we present a simple Lagrangian that generates the mixing matrices 
considered in the previous section. We find that the model naturally generates the 
BM mixing though the TBM mixing can also be obtained but with introducing 
different flavor symmetries. Our phenomenology will be done in the scenario in 
which the leptonic mixing is BM in the symmetric limit. 

We will use the seesaw mechanism to obtain the neutrino masses. Our model 
extends the SM by an additional Higgs doublet and two singlet scalars. The particle 
content of the model is given as 

• three left-handed lepton doublets Da^, where a denotes e, /i, and r, 

• three right-handed charged-lepton singlet a^, and 

• three right-handed neutrino singlets i/aR- 
In the scalar sector, we employ 



M 




(16) 
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• two Higgs doublets (pj with vacuum expectation values, v.e.v, (O|0°|O) = ^ 
and 

• two real singlet scalar fields ei and £2, with v.e.v's (0|e°|0) = Wk- 
The symmetries of the Lagrangian are introduced as 

Dei ^ei,, GR, l^eR ^eR, 

ei -ei, 62 62, 01 01, 02 02, 

■^2 • '^eR, gr, -DeL) ^ii ^2, (Change sign, and the rest of the fields remain same) 
Z2 '■ ^eRi ^fxR, i^tRi 02, (Change sign, and the rest of the fields remain same). 

(17) 

The most general Lagrangian consistent with the symmetries is 

Cy = [viD^^eR + y2 {D^^jiR + D^^tr) + {D^^tr + D^^hr)] 0i 

+ [yADej^Pe.R + Vb {D^^Vf^R + L'^i'^ril)] 02 + h.C, (18) 

Cm ^ \ [Mi/J^C- + Mpiy^jiC-\R + Mpu^j^C-^UrR] 

1 T r^-i f (aei + fees) {aei - be2)\ , . 

- i^V^crC [i^nR ^ + i^tR ^ 1 + h.c. (19) 

Here, 0^ = icr20* is the conjugate Higgs doublet and we have chosen to work in a 
basis where the Dirac mass matrix for the neutrinos is diagonal. We can simplify the 
Lagrangian in several ways. First, we can redefine aei — >■ ei and 6e2 — >■ 62. Second, 
to reduce the number of parameters we can impose an approximate symmetry of 
the Lagrangian. A SU{'i) symmetry where the right handed singlet fields and the 
left handed doublet fields transform as the SU{3) triplets leads to y^ = y^ = yo- 
The SU{3) symmetry is only satisfied by the Dirac mass term for the neutrinos 
and is broken by the other terms in the Lagrangian. Third, we will require the 
Lagrangian to be invariant under the transformation of the right-handed charged 
leptons {hr •(-)■ —trjCr — )■ —eR^cpi — 0i), with all other fields remaining un- 
changed. This symmetry requires 2/2 = Z/3 leading to vanishing ^ mass. The ^ mass 
is introduced later as a symmetry breaking term. Finally, we will set the Majorana 
mass terms M — Mp. We can then rewrite the Lagrangian as 

Cy = [yiDe^ep + 1/2 {D^^ixr + D^^tr) + y2 {D^^tr + D^^/ir)] 0i 

+ yo [DcL^eR + D^^V^R + Drj^UrR] 02 + h.c, (20) 
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- I^V^crC y^^R ^ + ^tR ^ 1 + h.c. (21) 

The most general scalar potential V that is invariant under X X is 
given by 

V = -/x?e? - + Aiet + A26^ + A;6?6^ 

+ (7iei|(/>l|^ + (726^1021^ + CT3e^|0l|^ + (14621021^ + ^2f/D(</'l, 02), (22) 

where y2ifD(0i, ^2) is the potential of the two Higgs doublets, 

+ ^^12 ('^I'^l - 4<^2)^ + A^,, ((0l0l)(402) - (0l</>2)(0^<?'l)) 

+ ((0l0l)(</'2</'2) + (0l02)(</'^</'l)) . (23) 

If we impose an additional symmetry to the above potential such as ei -H- £2, then 
the potential takes the form 

V = V(e? + 6^) + (e? + e^)X:a40. + A(e? + e^)^ 

i=l 

+ A'(e?-e9' + 1/2i/D(0i, 02). (24) 
We can parametrize the v.e.v's of the singlet scalars as follows 

(0|ei|0) = ^008 7 and (0 |e2| 0) = 'u;sin7. (25) 
Thus, the only term that depends on 7 is 

/(7) = AVcos^27. (26) 

By minimizing /(7), one gets 

cos 27 = 0. 



Thus 



(27) 



(0|6i|0) = (0|62|0) = ^. (28) 



By minimizing the above potential one can find the parameter w and the v.e.v's of 
the two Higgs doublets which are nonzero and different in the symmetric limit 



2ft' 



V. = ,/^, (29) 



where 



- 2A02/iVi + 2A'^^^/iV2 + 2A<^,2M^^2 + fJ'%(^i(^2 - - 2A'^^2/^^(^i + ^2), 

Pi = 4A( — A'l^i — 2A'^2l'^<?f'12 + ^<^l^<^12 + ^<^l^<ji2 + ^012-^02 + '^'(/>12(2'^'(/.21 + "^^l + ^^<i>12 

+ ^.^2)) - ^-i'12'^i - ^H'^i + 2A'^2i^i^2 + 2A<^,2aio-2 - A^.CTa - A0,2'^2 " ^12^^^ + ^2)^ 

Q;2 = 4A(A</,,2/i^2 + '^4>lf4>2 ~ -^(/>21^ll ~ ^<^12M|i + Kl2(f4>2 + /^ll ) ) ~ 2A<^i2/^^'^2 

- 2A0,/iV2 + 2A'^2^/iVi + 2A<^i2/^^^^i + /^Ji(^2(Ti - /ij2f^? - 2A'^^2/i^(a2 + ai), 

^2 = 4A(-A'|21 ~ 2A'^2l'^012 + ^<t>2^4>12 + -^02^01 + ^<t>12^4>l + •^'(/.12(2'^'(/.21 + "^02 + 4A<^i2 

+ '^4>l)) ~ ^^12'^2 - ^<l>i'^2 + 2A^2i'^20"l + 2A0i2<^2C^l - A<^2C7i - Xct>i2'^f - ^^012(^2 + CTif. 

(30) 

Also, the parameter w can simply be written as follows 

which shows that the v.e.v of the singlet scalars is independent of {vi, V2) when 

(7i = (72 = 0. 

The explicit form of the Yukawa matrix, Y^, and the Dirac neutrino mass matrix 
can be written from the Lagrangian (18) as follows 



'yi 

y2 y2 
y2 y2 



>^2i = ^ I y2 y2 I , (32) 



Md = diag(^. A, A), with A = y-^. (33) 

V2 

Also, the Majorana mass matrix can be obtained from Eq. (19) as follows 

M 

Mr= \ -v^ M I , (34) 





with Vu, — yw. Using the seesaw formula [18], the neutrino mass matrix is given as 

= -MIM^^Md. (35) 



Then M.^ has the structure 



where 



X 


G \ 


G 


X 





Z ) 








M2 - vl 



(36) 



X - G--^^ Z--^" (37) 

By diagonahzing Eq. 36, we obtain the neutrino masses as 



nil 

1712 



M + vJ 
'M-vJ 

\2 



= (38) 

Note that from the above equations one can estimate the scale of the v.e.v, V2, 
of the second Higgs doublet 02- As the absolute neutrino masses are in the eV scale, 
therefore, V2 has to be in the MeV scale if the see-saw scale (M) is in the TeV range. 
The mass relations satisfy the relation 

112 
— + — = —. 39 

mi 1712 1713 

Similar relations among the masses are discussed in Ref. [19]. We can use the above 
sum-rule to obtain an upper limit for the heaviest mass, ImQl 

< for the 

) I oi ^ ||miH-|m2|| 

normal hierarchy or \m2\ ^ for the inverted hierarchy. 



4 Symmetry Breaking 

The breaking of the flavor symmetries in the charged lepton and the neutrino sectors 
will cause deviation from the BM form, and we study these deviations in this section. 
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4.1 Charged Lepton Sector 

In the charged lepton sector we break the decoupled 2 — 3 symmetry by adding the 
following higher dimensional terms 

Oi = cy2D^,iiRct>^^, (40) 

and 

O2 = y' {De^fiR - D^^tr + D^^cr - D^^cr) 01^^- (41) 

The operator O2 breaks the decoupled 2 — 3 symmetry, x Z|, but is still 2 — 
3 symmetric. The operator Oi explicitly breaks the 2 — 3 symmetry, Zf^, and 
generates the muon mass. To generate explicit 2 — 3 breaking we have introduced 
the higher dimensional operator in the position of the muon field, 2-2 element, in the 
Yukawa matrix which is the most straightforward way to generate the muon mass. 
Introducing this operator in the 3-3 position generates the same numerical solutions 
for the correction angles. But introducing it in the 2-3 or 3-2 positions does not 
generate physical values for the mixing angles. Even introducing 2 — 3 symmetric 
terms in (2-2, 3-3) or (2-3, 3-2) generates either unphysical mixing angles or gives 
very large correction mixing angles that do not lead to successful phenomenology. 

In the presence of the higher dimensional terms the charged lepton Yukawa 
matrix has the following form 



^11 ^12 —^12 

= ( 1,2 IIt{1^+2ki) i/T I , (42) 

"^12 2^T 2^T 



with Ki = cvl/2A^ and /12 = y'vf/2\/2A^ after the Higgs field gets it's v.e.v. Three 
relations can be obtained among the matrix elements 



^12 — ^13^ 
yL _ y 
^23 ~ ^33^ 



L 

L 



F22 = (1 + 2«;0F23- (43) 

We can solve for the unitary matrix, [//, that diagonalizes Y^ in Eq. 42. We 
write, 

Ui = W^,R^,,R[,R[„ (44) 

where 



R12 



C121 = cos6li2i;si2i = sin6li2/, (45) 
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R 



13 



Cm = 



1 
-sisiC^* Cm 

cos9m]Sm = sin6'i3;, 



(46) 




-^23 - 

C23« = COS 6^23; ; S23i = sin 6123/ • (47) 
The Yukawa matrix, Y^, can be written as 

= c/,y,^c//, (48) 

with 

(le 0\ 

yj^ = I, . (49) 

\0 Ir) 

Applying the relations in Eq. 43 to the matrix elements in Eq. 48 using Eq. 44, 

one can solve for the corrections of the mixing angles. Two ways can be used to find 
the angles, analytically or numerically. Solving for the mixing angles analytically, 
see details in appendix A, can determine the size of the Yukawa matrix parameters 
in Eq. 42 



rrir ' 

Ki — ZZf^, 



12 



It « {lr-W-\{zz,n (50) 
with z ^ 2. It is interesting to note that 

K^^^ZZ^ (51) 

which fixes A ~ TeV. We assume that the charged lepton corrections are "CKM- 
like", i.e. 

sin^i2/~A, sin 6*23/ ~ AA^, sin ^13; ?s SA^, (52) 

where A and B are real and of order one and A is the sine of the Cabibbo angle, 
A = sin^^c — 0.227. Wc present numerical solutions to the relations in Eq. 43 for 
various z values that produce the pattern in Eq. 52. In our calculations we assume 

5 — TT, 
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• For z = 2.0: S121 ^ ±0.34, S131 ^ ±0.0011, S231 ^ -0.059, 

• For z = 2.06: sui ^ ±0.3, S13/ ^ ±0.001, S231 ^ -0.061, 

• For ^ = 2.2: Sia ^ ±0.2, Si3, f« ±0.00075, 533; ~ -0.065, 

We expand the angles in Eq. 1 as 

r 1 1 

«i3 = ^, si2 = ;^(l + s), 523 = ^(1 + a), (53) 

where the three real parameters r, s, a describe the deviations of the reactor, solar, 
and atmospheric angles from their bimaximal values. We use global fits of the 
conventional mixing parameters (s, a) [20] that can be translated into 3a ranges 
and the mixing parameter r with 2.5cr significance (90% C.L.) [5] 

0.12 < r < 0.39, -0.29 < s < -0.14, -0.15 <a< 0.16. (54) 

To first order in r,s,a, the lepton mixing matrix can be written as. 



U 



-i(l + s-a + ^e^^) l(l-s-a-^e^^) ^(1 + a) |, (55) 

which is similar to the parametrization in Ref. [21] with the TBM mixing. We have 
assumed that 5 — tt where the present data prefers a negative value for s [21] and 

r is positive, in our discussion we do not consider CP violation. Now, we can write 
the parameters (r, s, a) in terms of the elements of the mixing matrix 

s = -I + V2U12, 

r = V2{l + s-a + 2U2i), 

a = -I + V2U23. (56) 
From the details in appendix A, one obtains 

S ~ -^{Sl2l + Si3l), 

r ^ S121-S131, 

a w -S231- (57) 
From the above equations one can get the deviation parameters as follows 

• For z = 2.0: s f» -0.24, r ^ 0.34, a ^ 0.059, 

• For z = 2.06: s f« -0.21, r ^ 0.30, a 0.061, 

• For ^ = 2.2: s ft; -0.14, r ^ 0.20, a 0.065. 

The above results demonstrate that the contributions from the charged lepton sector 
can accommodate the T2K data of 613 as well as the other mixing angles. 
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4.2 Neutrino Sector 



In this section we consider deviations of the BM mixing from the neutrino sector. 
We maintain the invariance of the Majorana Lagrangian under the symmetry group 
in Eq. 17 and generate the deviation from the BM matrix by breaking the ei ■H- 
62 symmetry in Eq. 24 by introducing the most general dimension four symmetry 
breaking terms in the potential 

{el-el)i2^m + e{el-el){el + el)^ (58) 

i=l 

We require that all terms in the symmetry breaking potential are of the same size 
which results in p ~ '^^i where v is the electroweak v.e.v with v'^ — vf + and ui 
is the scale of the v.e.v's of the singlet scalars. Thus, the potential is 

V = {el + 6^) + (6? + 6^) a,<P% + A (6? + 6^) V A' {el - 6^) ' 

i=l 

2 

+ (ei - £2) E '^^^^'^^ + q{^1- 4) {4 + el) + V2hd{<Pi, <p2). (59) 

i=l 

Now, parameterizing the v.e.v's as in Eq. 25 and minimizing the potential leads to 

Qw'^ + (cril^ip + a2\v2\'^) 



cos 27 = 



2A'w2 



^2 ^ 2AX//^ - (ail^ip + ^21^^21')) + q{<^'i\vi\'' + (^2^21') ^gQ^ 

Keeping in mind the size of the various co-efficients in the symmetry breaking po- 
tential discussed above, we find that cos 27 ^ up to corrections of order We 
assume that w is in the TeV scale and with v in the EW scale the symmetry breaking 
corrections are of the right size to explain the experimental numbers. 

We shift the v.e.v's of the two singlet scalars {wi ^ W2) up to the first order of 
the symmetry breaking parameter. Then, the Majorana neutrino mass matrix in 
Eq. 34 takes the form 

Mr = I -v^p M I , (61) 
-t;,™, M 



where 



Vwp = -j={Wi+W2), 



-l=(w,-W2). (62) 



^The most general symmetry breaking terms can be expressed in terms of the form in Eq. 58 
and symmetry conserving terms that can be absorbed in the symmetric potential. 
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We write the v.e.v's of the singlet scalars after symmetry breaking as 



w + pi 



V2 ' 

W + P2 .^oN 

W2 = ^ , (63) 



where pi and p2 are small quantities and 



WT 

Pi = -p2 = —■ (64) 



Up to the first order of the symmetry breaking parameter r, 



^ = ' ^^^^ 



one gets 



T 

Vwn = 2^"'- 

It turns out that breaking the ei •<-> £2 symmetry to generate different v.e.v's for the 

singlet scalars is not sufficient to break the almost degeneracy of (mi, 1712) to satisfy 
the squared mass difference measurements. Therefore, we introduce an additional 
term in the Lagrangian which is consistent with the symmetries of the Lagrangian, 

Ml Kr<^-\r + '^JrC-'u^r] . (67) 

Thus 

(68) 





/ M 




-\v, 


Mr= I 


-Viu 


M' 












M' 



where M' = M + Mi. 

The neutrino mass matrix in Eq. 36 changes to be 




= \ G' Y' W \ , (69) 

where 



X' = 
Y' = 



AA^M' 



A^{AMM' -vIt 



M'{AMM' -vI^A + t"^))' 
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Z' 

a 
p' 

w 



M'(4MM' -^2(4 + ^-2))' 



4MM' -172(4 + r2)' 



4MM' -'y2(4 + T2)' 
M'{AMM' - vl{A + t2)) ' 



(70) 



By diagonalizing Eq. 69, one gets the mass eigenvalues 

2^2 ({M + M') - yM2 - 2MM' + M'2 + 'f;2 (4 + t2)) 

4MM' -'i;2(4 + t2) 
2^2 J|(M + + ^M'^ - 2MM' + M'2 + 't;2 (4 + r^j^ 



mi 



m2 



4MM'-i;2(4 + t2) 

".3 = (71) 

Now, we can diagonalize the mass matrix in Eq. 69 using the unitary matrix 
= W^2R23Ri2 with, 

ij 

Ci2u = cos 9i2u; Si2u = Siu 9i2u, 

^ _ /I M 

-^23 ~ I '^231^ -523,/ I , 
\0 — S23^ C23;/ / 

C23I. = COs6'23iy;S23iy = Sin6'23i/. (72) 

The mass matrix elements in Eq. 70 satisfy the two relations 

X\Z'-Y') = P'^-G'\ 
G'P\Z'-Y') = W\P'^-G'''). (73) 

By applying the above relations to the matrix elements of 

= U^MtUl (74) 
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one can obtain the mixing angles 



'2mi(m2 — rris) 
m2(mi - ms) ' 



, -mim2 + 2mim^ - 7712171^ 

si2u = \ ^ ■ 75 

im^yrai — 777.2) 

Eventually, we obtain the elements of the lepton mixing matrix Upmns = UjUi, 
with = 1^23-^23-^13-^12 Ui, = Wi2R2zRi2- The deviation parameters (s, r, a) 
can be obtained from Eq. 56 as follows 

S ~ J=\Si2l + Si2,i) + Si2i„ 

v2 

r ~ S121 - Si3i - S22,v, 

a ~ + ;^'523i.- (76) 



5 Numerical Results 

Prom the neutrino mass matrix (16), one observes that in the degenerate case, when 
7773, a c, d Q which means that the neutrino mass matrix is already 
diagonalized as M.^, ~ diag (a, a, a). That means the lepton mixing matrix does not 
include a contribution from the neutrino sector, and the resultant leptonic mixing 
is inconsistent with the experimental data. Thus, in the symmetric limit our model 
excludes the case of the degenerate neutrino masses. Even, after symmetry breaking, 
the degenerate case in Eq. 71 leads to vanishing the v.e.v's of the singlet scalar fields 
which does not lead to successful phenomenology. 

The numerics goes as following; we choose masses (7771, 7772, 7773) which satisfy the 
experimental values of the squared mass differences 

A777^i = ml-m\^ (7.59 ± 0.20) x IQ-^eV'^, 

A777^2 = |777^ - 777^1 = (2.43 ± 0.13) X lO^^e^^ (77) 

We substitute those mass values in (r, s, a) in Eq. 56, using (si2i/) ^2Zv) given in 
Eq. 75 and (si2«, S23«, Si-n) in sec. (4.1). If the results satisfy the experimental 
constraints in Eq. 54, we plot the possible values of the absohitc masses and the 
mixing angles. By using Eq. 71, we calculate values for the Lagrangian parameters 
[vw, A, M, M') which generate the values of the absolute masses obtained from 
the graphs. Prom the graphs, one find that {vy,, M, M') are obtained in the TeV 
scale and A in the MeV range. 

Three mass-dependent neutrino observables are probed in different types of ex- 
periments. The sum of absolute neutrino masses rricosm = ^Trii is probed in cosmol- 
ogy, the kinetic electron neutrino mass in beta decay (Mg) is probed in direct search 
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for neutrino masses, and the effective mass (Afge) is probed in neutrinoless double 
beta decay experiments with the decay rate for the process T oc M^g. In terms of 
the "bare" physical parameters rrii and Uai, the observables are given by [19] 

Urui = \mi\ + \m2\ + jmal, 

= ||mi||[/elP + 1 1 f/e2 1 ^-^^ + | 1 1 f/e3 1 | , 

= V|mi|2|[/eiP + |m2|2|C/e2|2 + ImaPlC/gaP. (78) 

In our analysis we ignore the Majorana phases (0i, ^2) and plot Mg versus Em^ 
and versus mught, where mught is the lightest neutrino mass. 

In Figs. (1, 2, 3) we assume specific values of z with the corresponding correction 
mixing angles (sia;, S121, S23/) and plot the absolute masses and the mixing angles 
which satisfy the neutrino mixing constraints. By choosing a value for the symmetry 
breaking term r, we plot the parameters (u^, A, M, M') that satisfy the squared 
mass difference measurements. This model supports the normal mass hierarchy as 
shown in the graphs with the scale of the neutrino masses in the few meV to ~ 50 
meV range. The results agree with the recent T2K data which find a relatively large 
^13. The graphs show that the see-saw scale (M, M') are in the TeV range, and 
the second Higgs that couples to the right-handed neutrinos has v.e.v ^2, included 
in A, in the MeV scale. Also, they indicate that the v.e.v of the singlet scalar fields 
Vw is in the TeV scale. The graphs show that Em^ ~ 0.06 eV and Mee < Mg and 
Mee < 0.35 eV [22]. Various other mechanisms to generate the neutrino masses with 
TeV scale new physics are mentioned in Ref . [23] . 

6 Conclusion 

In this paper we presented a model for leptonic mixing which accommodates the 
sizable neutrino mixing angle ^13, recently measured by the T2K and MINOS exper- 
iments. We worked in a basis where the charged lepton mass matrix is not diagonal 
and proposed an explicit structure for the charged lepton mass matrix which is 2-3 
symmetric except for a single breaking of this symmetry by the muon mass. We 
identified a fiavor symmetric limit for the mass matrices where the first generation 
is decoupled from the other two in the charged lepton sector while in the neutrino 
sector the third generation is decoupled from the first two generations. The lep- 
tonic mixing in the symmetric limit was shown to have, among other structures, the 
bi-maximal (BM) and the tri-bimaximal (TBM) mixing. 

A model that extended the SM by three right handed neutrinos, an extra Higgs 
doublet, and two singlet scalars was introduced to generate the leptonic mixing. 
In the symmetric limit the model had two Z2 symmetries in addition to the n — r 
symmetry and the BM leptonic mixing was obtained when the two singlet scalars 
got equal v.e.v's. 
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Symmetry breaking effects were included in the charged lepton sector via higher 
dimensional operators that generated a — r symmetric mass matrix except for 
a single breaking due to the finite muon mass. In the neutrino sector, symmetry 
breaking was included via shghtly different v.e.v's for the two singlet scalars. To 
explain the Am^ data two different Majorana mass terms, one for and one for 
and Vri was used keeping in mind that the fJi — T symmetry fixes the Majorana mass 
terms for the and to be the same. 

A fit to the experimental measurements showed that our model predicted normal 
hierarchy for the neutrino masses with the masses being in the few meV to ~ 50 
meV range. The Majorana mass terms as well as the v.e.v's of the singlet scalar 
fields were predicted to be in the TeV scale and consequently the v.e.v of the second 
Higgs doublet was shown to be in the MeV range. We calculated predictions for the 
mass-dependent observables (Smj), (M^) and (Mge). We found that Em^ 0.06 
eV, Mee < M^, and Mee < 0.35 eV. 
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Appendices 



A Charged lepton sector 

We analytically calculate the deviation of the leptonic mixing from the symmetric 
limit due to corrections from the charged lepton sector. We, here, are going to 
determine the sizes for the Yukawa matrix elements in in Eq. 42. We first 
consider the breaking of the 2-3 symmetry in the charged lepton sector via the 
introduction of a higher dimensional operator that generates the muon mass 

Oi = cy2D^^fXR(f)i^^jp-. (79) 

Thus, we consider the Yukawa matrix, 

fie \ 

= |/T(l + 2fi:0 . (80) 

\0 ¥t \It) 

The structure above breaks the 2 — 3 symmetry because of the correction to the 22 
element. Note that wc do not break the 2 — 3 symmetry in the 23 element so that 
the Yukawa matrix remains symmetric. The matrix is now diagonalized by the 
unitary matrix, Ui = W22R23- Applying the relation (^23) 23 = (^23)33 Eq. 80 to 
= UiYiagUl leads to 



t 



231 



2 



^ ■ " ^ (81) 



where t23i = tan ^23; and we have chosen the solution that leads to small angle ^23^ 
and to small flavor symmetry breaking. Keeping terms up to first order in we get 

^23/ ~ -^M- (82) 

We further obtain for ni and It in Eq. 80, 

Ki ^ - tan 26l23i ~ 2^;^, 

h = (Z^-y cos 2^23/- (83) 

Comparing the above equation with Eq. 79, the size of the higher dimensional op- 
erator can be estimated as 

^2 - 2.,. (84) 
Since Vi ~ 250 GeV therefore the scale of A is in the TeV range. 
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To obtain a realistic charged lepton matrix, we take into account the mixing 
involving the second and the third generations in the full Yukawa matrix 

O2 = y' {D.^fiR - De.TR + D^^en - Dr^en) (85) 
Thus, the full mixing matrix will be give by 

(^11 ^12 —^12 \ 

h2 |Zt(1 + 2«:0 . (86) 

We will assume that the Yukawa matrix is now diagonalized by the unitary 
matrix Ui given by 

= ^^23-^23-^13-^12- (87) 

Prom the Yukawa matrix (86), one can find the two relations 

^12 = ~^13) 

Y22 = ^(F23 + 1^33)(l + 2«:;)- 



Applying the above two relations to 

■ diag I 



= UiY^^u; (89) 



using Eq. 87, one can obtain the solutions 

S12I ±C23i 



1^ - 2ki + (-3 + SZf, - 2ki)c231S231 + 2Zi^KiC23lS23l 



- + 2k«)c23«-S23Z 



, ^ Z^- 2ki + (-3 + 3Zf, - 2ki)c231S231 + ^Z^KiC23lS23l .„„^ 

sm ~ ±^/^^ci2iS23i\ — ^^^^ ^ - 90 

*' Z,_,{2> + 2ki)C231S231 

By comparing Eqs. (86, 89), one can get the matrix element I12 after substituting 
Eqs. (82, 83, 90) up to the first order in 



^12«yf (/e-g. (91) 

The leptonic mixing matrix is now given by 

UpMNS ^U}U,, (92) 
where Ue = W^23^23^L^'i2 and = W^^- 
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Figure 1: Scatter plots for z = 2.0 with S12/ ~ —0.34, Sm ~ —0.0011, and S23/ 
—0.059. In the neutrino sector, we assume that r = 0.1. (meV = 10^^ eV) 
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Figure 2: Scatter plots for z = 2.06 wit4P Si2i ~ —0.3, S13; ^ —0.001, and 523; 
—0.061. In the neutrino sector, we assume that r = 0.05. (meV = 10~^ eV) 
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Figure 3: Scatter plots for z = 2.2 with25'i2z ~ —0.2, Si^i ~ —0.00075, and S23/ 
—0.065. In the neutrino sector, we assume that r = 0.1. (meV = 10^^ eV) 



